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Abstract 

We combine Maxwell's equations with Eulers's equation, related to a 
velocity field of an immaterial fluid, where the density of mass is replaced 
by a charge density. We come out with a differential system able to 
describe a relevant quantity of electromagnetic phenomena, ranging from 
classical dipole waves to solitary wave-packets with compact support. The 
clue is the construction of an energy tensor summing up both the elec- 
tromagnetic stress and a suitable mass tensor. With this right-hand side, 
explicit solutions of the full Einstein's equation are computed for a wide 
class of wave phenomena. Since our electromagnetic waves may behave 
and interact exactly as a material fluid, they can create vortex structures. 
We then explicitly analyze some vortex ring configurations and examine 
the possibility to build a model for the electron. 
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1 Introduction 

Since the advent of the theory of electromagnetic fields, more than a century 
ago, waves have been described as a kind of energy flow, governed by suitable 
transport equations in vector form, namely Maxwell's equations. The original 
publications (see [19] and [20] ^ . in a way that nowadays may seems naif, ac- 
tually refer to electric and magnetic phenomena using mechanical terms, such 
as tension, stress, pressure and vortices, in association with some 'imaginary 
substance'. 
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It turns out that, in void, the electric and magnetic fields (E and B, respec- 
tively) are transversally oriented with respect to the direction of propagation, 
and their envelope produces a sequence of wave-fronts. This is in agreement 
with the fact that the energy develops according to the evolution of the vector 
product E x B, otherwise known as Poynting vector. 

On the other hand, the dynamical behavior of a compressible non viscous 
fluid is well described by Euler's equation, where, in principle, the velocity vec- 
tor field (denoted by V) might not be necessarily related to a real material 
fluid. In particular, one could replace the mass density by a sort of charge den- 
sity. Therefore, the temptation to describe electromagnetic and velocity fields, 
through a combination of the respective modelling equations, is well motivated. 

We are going to present a system of equations in the three independent 
vector unknowns: (E, B, V). In pure void, the electric and magnetic fields follow 
the Faraday's law together with the Ampere's law, where a current, flowing at 
velocity V, is supposed to be naturally associated with the wave. In order to 
close the system, the third relation is the Euler's equation for V, containing an 
added forcing term E + V x B, perfectly analogous to that characterizing the 
Lorentz's law. In this way, the three entities (E,B,V) turn out to be strictly 
entangled. 

Despite the appearance, the new model allows for a very large space of so- 
lutions. Moreover, it displays numerous conservation and invariance properties, 
all deducible from a standard analysis. An interesting invariant subspace of so- 
lutions is the one where the third equation is reduced toE + VxB = 0, which 
means that no acceleration is acting on the wave, and the corresponding 'flow' 
is somehow laminar. For this circumstance, the solutions, called free-waves, 
perfectly follow the laws of geometrical optics, ruled by the eikonal equation. 
Together with other known solutions, free-waves also include solitary electro- 
magnetic waves with compact support almost of any shape, intensity, frequency 
and polarization. Such a result is clearly important, since it reopens the path 
to a serious discussion, based on a deterministic analysis, on crucial issues as 
the structure of photons, the duality wave-particle and the quantum properties 
of matter. 

Far more complicated solutions (not of the free- wave type) are however possi- 
ble. Since our electromagnetic radiations actually behave as a fluid, they can be 
constrained to evolve in bounded regions of space, similar for instance to vortex 
rings. According to the model equations, rotating photons in a vortex structure 
may carry a charge and deform, via Einstein's equation, the local geometry of 
the space-time in order to create a gravitational environment assimilable to the 
presence of mass. The same metric space is responsible for the stability of such 
a wave, obliged in this way to develop along self-created geodesies. This leaves 
us with the conjecture that some stable elementary particles (as for instance 
the electron) could be made by rotating photons, an idea that has been put 
forward by many authors in the past, although with not too much recognition, 
basically due to the lack of a sufficient theoretical description of electromagnetic 
phenomena, able to go beyond the classical linear Maxwellian approach. With 



2 



the help of the new material collected here, we examine this aspect in the last 
section. 

This paper is a review and a reorganization of the results presented in [11] , 
and contains additional new material. The presentation follows an inverse path, 
and the equations are first derived in covariant form and successively deduced 
in the flat space. The general relativity framework is however the key point to 
understand the structure of the solutions, especially in the case of constrained 
waves. The entire theory comes out as a consequence of a wise choice of the 
right-hand side tensor in the formulation of Einstein's equation. This combines 
in the proper manner the electromagnetic stress tensor and a mass type tensor, 
where the density is referred to the electric field. In the subspace of free- waves, 
by generalizing the results given in [11] , we are able to write down explicit 
general solutions of the Einstein's equation. As a final step, we will derive a 
notion of mass density, that is going to be zero in the average for free-waves 
(and photons) and positive for electrons. 



2 The right-hand side tensor 

Let us fix the notation by introducing the Einstein's equation: 

Rap — \g a pR = — xT a p (1) 

where R a p is the Ricci tensor, R the scalar curvature (the trace of R a p) and x 
a positive adimensional constant. The signature of the metric space g a p will be 
(+,—,—,—). The given constant \i that will be estimated at the end of section 
6, has nothing in common with the gravitational constant G. 

As anticipated in the introduction, the main ingredient of our theory, is the 
choice of the right-hand side tensor. We set: 

T al3 = liUctfl + Map) (2) 

where c is the speed of light and fj, is a constant, dimensionally equivalent to a 
charge divided by a mass (see (|62jl ). In ([2]) we find the usual electromagnetic 
stress tensor: 

U aP = -9 7S F a ^Fp S + y a pF j5 F^ (3) 

with 

F aP = V a A{i-VpA a (4) 

where A a is the electromagnetic 4- vector potential (energy/charge) and V Q de- 
notes the covariant derivative. Although it is not strictly requested, A a satisfies 
the Lorenz gauge relation: 

V a A a = (5) 
In @ we also find a kind of mass tensor: 

M a p = pV a Vp - g a pp (6) 
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where V a is velocity 4-vector, p plays the role of a density and p is similar to a 
pressure. Dimensionally, p and p are taken in order to conform with the setting 
in @, therefore they will not correspond to the usual entities of classical fluid 
dynamics. Moreover, differently from what it is generally supposed (see for 
instance [TUj, P-91), no direct relations between p and p are a priori requested. 

This way of building T a p characterizes the whole theory. We remark (and 
this is an important key point) that U a p in ^ appears with the opposite sign. 
The reason for this choice will be clear as we proceed in the exposition. 

Concerning the vector field V a , we shall make the following assumption: 

g a0 V a V p = (7) 

that is equivalent to an eikonal equation for suitable propagating wave-fronts. 

A first important relation is obtained by evaluating the trace of (TT]). This 
gives: 

P = -^—R (8) 

In computing (JHJ) we considered that both the traces of U a p and pV a Vp are zero. 
In the first case this is a known property of the electromagnetic stress tensor. 
In the second case it is a straightforward consequence of relation ([7]) . Equation 
ijHJ) relates curvature and pressure. Wave-fronts travelling unperturbed, with- 
out developing interstitial pressure, actually move in a space with zero scalar 
curvature, and vice versa (see section 4). 

Let us then derive from {T]) a set of covariant equations. In order to have 
compatibility with the left-hand side in (TIJ, the most important relation to be 
satisfied is: 

VpT a P = (9) 

yielding the set of modelling equations. Due to the special expression of T a p, 
the energy conservation property ^ mixes up electromagnetic and mechanical 
terms in a very natural way. Thus, we start by observing that (see for instance 
[10] or [22]): 

V p U ap = g^F^(V/sF^) (10) 
V M afj = P G a + V f3 (pV l3 )V a - g^Vpp with G" = 1/ 7 V 7 V Q (11) 
Thus, according to ([H]), by summing up the above expressions, one must have: 

- ng^F^CVpFtP) + P G a + V„(pV^)V a - g^Vpp = (12) 

By adding and subtracting the term (pp/c)g aj F 7 sV s — (pp/c)F aS Vg, equation 
IT2l) becomes: 



+ V (pV )V a - g aP V f3 p = (13) 
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Hence, in order to obtain (|13[) . a sufficient condition is to impose the following 
set of equations: 

c\7pF s/;l = - pV s (14) 

p(G a + ^F aS V s ) = V a P (15) 

that, being written in covariant form, are, in principle, admissible in any geomet- 
ric environment, but, hopefully, in the metric spaces compatible with relation 
([7]). Equation fT4")) is the Ampere law and (Ti"5"j) is the Euler's equation for a 
certain velocity field (not corresponding to a material fluid). In truth, in order 
to get (fl"3|) . one should also impose the following continuity equation: 

V fi {pV fi ) = (16) 

but this is a consequence of equation (fl4j) . In fact, by taking the 4-divergence of 
(fT4|) . one has: cV^V^F" 5 ' 3 ) = = Vs(pV s ), resulting from the anti-symmetry 
of the tensor F Sf3 '. 

There is an interesting compatibility relation to be remarked. First, thanks 
to (jTj) , one recovers that: 

G a V a = \{G a V a + G a V a ) = i(^V 7 )(y Q t/ Q ) = (17) 

Then, due to the anti-symmetry of F 5 @ , (fTSf implies: 

V a V a p = (18) 

which is a sort of conservation property for p (a kind of Bernoulli's principle, if 
one considers that p = p + \ pV a V a ) . 

The equations (fT4"|) and (jT5]l . together with the equation of state ((5J), are the 
foundations of our theory. In principle, one could obtain §§§ without necessarily 
enforcing (fT4"|) and (|15p . However, we believe that all the most relevant phe- 
nomena can be described through the two differential equations written above. 



3 The model equations in Minkowski space 

To better understand what is happening we can specialize equations (fT4")) and 
(fT5|) in the case of a flat metric space g a p = diagjl, —1,-1,-1}, in the Carte- 
sian coordinate system (ct, x-y, x%, x^). 

We start from the potential (A ,Ai,A 2 ,A s ) = (— $,A) and through (g]), 
considering that V a = {\-§i, V), we obtain the electromagnetic tensor (see [TU] . 
section 24, or pj, section 5.2): 



/ 





-E x 


-E 2 


-E 3 
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-cB 2 






E 2 
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(19) 
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where E = (Ei, E 2 , E 3 ) is the standard electric field and B = {B\, B 2 , B 3 ) is 
the magnetic field. The above construction is true by virtue of the setting (see 
Q): 

B = - cmiA E = V$ (20) 



c 

In particular, we have: 

dB 

~dt 



dt 



curlE divB = (21) 



as a direct consequence of (f20| . The contravariant tensor F a @ is obtained from 
F a /3 by replacing E by — E. 

Without loss of generality, we can set V° = — c, so that, from (| 14[) with 
S = 0, we deduce: p = V /3 F ' 3 = divE. Denoting by V = (Vi,V 2 ,V 3 ) the 
velocity field, by one obtains V a = (-c,V) and V a = (-c,-V). Finally 
we get: 

<9E 

— = c 2 curlB - pV (22) 

% = ME-V) (24) 

where j^V = ^V+(V- V)V is the substantial derivative, so that (|23|) ends up 
to be the Euler's equation. Here equation (|22|) comes from (fT4|) for 6 = 1,2,3, 
equation |23|) comes from (I15|) for a = 1,2,3, and equation (|24|) comes from 
(fT5|) for a = 0. In addition, one has the continuity equation: 

|! = -div(pV) (25) 



P \~Dt + ME + VxB)) = -Vj> (23) 



that is easily derived either from (|16|) or by taking the divergence of (|22j) . 

The equations (jSJ) and |7]) are not valid in this context, since they are strongly 
related to the fact that g a p is solution of (flj , which is not true in the simplified 
case we are examining. The solutions to (|22l) - (|2"4")l are however important to 
set up the tensor F a p in Q, which may be computed using partial derivatives 
in place of covariant derivatives. This enables us to build the tensor |J5J| for a 
generic g a p to be used as a right-hand side in (TTJ). 

An even more simplified version is obtained by requiring DY /Dt = and 
p = in (j23|) : 

<9E 

— = c 2 curlB - pV (26) 

— = -curlE (27) 
divB = (28) 
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E + V x B = (29) 

where p = divE, and V is a velocity vector field satisfying |V| = c. Note that 
(J7J) is now true. The held V is oriented as the vector held E x B. Now © is also 
correct, since the property R = is compatible with the flatness of the space. 
Note that relation ([29")) is certainly satisfied for all electromagnetic waves where 
E is orthogonal to B and |E| = |cB|. Moreover, from (J29J we get E • V = 0, 
that is in agreement with ([24]) . since p is identically zero. Take into account, 
however, that the flat metric does not satisfy flj, since in this case one has 
R a /3 = 0, in contrast to the fact that T a p ^ 0. 

This set of equations produces a closed subset of solutions called free- waves, 
that contains all the electromagnetic phenomena in vacuum travelling undis- 
turbed according to the rules of geometrical optics. Indeed, if V = V\P is a 
gradient, then the relation |V| = |V^| = c is the eikonal equation. This sub- 
space includes solitary waves with compact support, perfect spherical waves, 
and many other solution non obtainable with the classical Maxwell's setting. In 
this case, equation (|26p comes from finding the stationary points of the standard 
Lagrangian of the electromagnetism, under the condition A a V a = on the 
potentials (see [14 ). Moreover, (|26p is invariant under Lorentz transformations 
(see [TT], section 2.6). 

The subset of free-waves is characterized in covariant form by the following 
expression, deduced from fTS]) for p = and G a = 0: 

F aS V s = (30) 



that brings to ([29]) and E ■ V = 0. It can also be proven that both ([25]) and J30J 
are invariant under Lorentz transformations. 

The triplet (E,B,V) turns out to be right-handed (in the right-handed 
reference frame (x,y,z)). There is however an alternative path, that brings to 
left-handed triplets. It is enough to replace ([20]) by: 

i i ft A 

B = - curlA E = V$ (31) 

c c at 

and rewrite all the equations with — B in place of B. We get in this way a 
completely new set of electromagnetic waves that have a specular image with 
respect to the classical one. 



4 Explicit solutions 

In the case of free- waves explicit expressions of the metric tensor g a p are avail- 
able. The results we are going to show generalize those given in [IT] . 

Let us assume the following form for the potentials: 

A a = (-/(i,yMO, 0, 0, f(x,y) T (0) (32) 
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Here we are in Cartesian coordinates (xq, X\, X2, X3) = (ct,x,y,z) and we set 
£ = z — ct. The two functions / and r are arbitrary. 

With this setting, we obtain the electromagnetic fields: 



E = (~f X T, - fyT, 0) 



B 



\fyT, ~f x T, 0) 



(33) 



where f x = df/dx and f y = df/dy. Thus, we are in presence of fronts parallel 
to the plane (x, y) and evolving in the direction of the z axis at speed c. Plane 
waves are obtained when / is of the form ax+by, for some constants a and b. The 
case when / has compact support is quite interesting, since it give rise to self- 
contained solitary waves, travelling straightly at the speed of light (photons). 
These solutions do not belong to the Maxwellian theory (see also [T3] and [Tl]). 
By defining V a = (Vb,V) = (— c, 0, 0,c), the entire set of equations is 
satisfied in the flat space. Thus, we can also write £ = z — V^t. We observe that 
divB is automatically zero, while we do not expect divE = —(f xx + fyy) T to 
be zero. 

For any arbitrary (nonzero) function s of the variable x + y, we look for a 
metric tensor having the following structure: 



5oo 



1 



.911 



2„2 



922 



912 = ff21 



333 



-fxfyCT 



1 



(34) 



being the other entries equal to zero. Here a is a function of the variable £. 
Note that, in this situation and in accordance to ([7]), we have: 



V a V a 



and V a A a 







(35) 



with V a = (-c,0,0,-c). 

The determinant of the metric tensor is: g 
reason, it is necessary that s ^ 0. We have: 



-[*(/* - fy)°} 2 - For this 



F, 



a/3 



( f X fy \ 

-f x f x 



-fy 



fy 



F 



a/3 _ 



( 



-1 



V ~f X -fy J 



1 \ 



10 0+1 

-10 0-1 

V -1 +1 j 



(36) 



where T = ±ts/<j, and the sign is set up according to that of s(f x — f y )cr. 
By explicit computation, one can get the Ricci tensor, where the nonzero 



entries are: 



Roo = R33 = —R30 = —R03 — — o" ja 



(37) 



In addition, one gets R — 0. Amazingly, the functions s and / disappear. For 
completeness, we also report the nonzero Christoffcl symbols, for s 2 constantly 
equal to 1: 

p0 -n3 f2„J -nO r 3 

1 11 = 1 11 = -Jx aa 1 22 - 1 22 - -fy acT 

r?2 = r^i = r^ 2 = = —f x fy<jcr 
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Tii = —^11 = fxx/ifx — fy) r3, 2 = ~ T 22 = fyy/(fx — fy) 

-^12 = T 2 i = — r^ 2 = — T 21 = fxy/{fx — fy) 
1 T^l _ r<l _ r>l T^2 r 2 r2 p2 _ fx a 



_"pl T^l T^l T^l T^Z T^Z T^Z T~V 

_i 01 ^ ~ 1 10 — 1 13 — 1 31 — — 1 31 — — 1 13 — 1 01 — 1 10 



- fy) 



T^l T"U T~ll T^l T^2 T^2 T^2 T^2 f / O Q \ 

— 1 02 ^ ~ 1 20 — 1 23 — 1 32 ~ — 1 32 ~ ~~ 1 23 — 1 02 — 1 20 — / r c \ \°°) 

V{Jx - fy) 

We should avoid the points (x, y) such that f x = f y , but this is in general a set 
of measure equal to zero. The coefficients of the metric tensor remain however 
smooth everywhere, thus the fact that the determinant g can be zero somewhere 
is not a crucial problem. 

On the other hand, based on the proposed metric, we can evaluate the 
electromagnetic stress tensor ([3]), where the nonzero entries are: 

U o = C/33 = -Uao = -U 03 = r 2 /a 2 (39) 
In the new metric space, one also discovers that: 

P = V,^ = ~j=£-W=9 F^) * + «T) = Q (40) 

V^9 dx P \F I 9 \ dx dy J 

where we used that s is a function of x + y and t/ct does not depend on x and 
y. Thus, we can get rid of the mass tensor in |J5J|. We can also check that (|30[) 
is true, ensuring that the concept of free- wave is, in some sense, invariant. As a 
further confirmation, one can also check that (see (fTTj) '): 



r)V a 

G a = v 13 —- + r™ s v~<v s = (41) 

OX/3 



because V a = (— c, 0, 0, — c) is constant and the Christoffel symbols Tq , Tq 3 , T 30 
and r^g are zero. The above properties are compatible with the choice p = 0. 

In the end, putting together (|3"T|) and (f3"9"]l . we find out that g a p is solution of 
the Einstein's equation (JT|) , if the following differential equation in the variable 
£ is verified: 

2 

II 2 (Ar,\ 

-ao = —t (42) 
c 

For example, when r(£) = sinw£, equation (|42p admits the simple solution given 
by: 

= ^ sincC (43) 

CTCJ 

Therefore, g a p is associated with a monochromatic gravitational wave strictly 
tight to the given electromagnetic wave, and having the same frequency and 
support. The intensity of a is inversely proportional to the frequency. This 
somehow agrees with the observation that gravitational phenomena are more 
relevant when one is dealing with very low frequencies. If we could ride a wave, 
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we would not be able to 'see' in the transversal direction, since, in the modified 
metric, the tensors R a p and U a p (see ([37|) and ([39]) . respectively) do not depend 
on x and y. This justifies why in this case we have p — (independently of 
the value divE originally attributed in the flat space). This means that, in its 
geometry, a photon looks like an unbounded plane wave. Note that here we 
solved the Einstein's equation in full form, and not its linearized version, as it 
is usually done in other contexts. As pointed out in [1] this is a decisive and 
necessary improvement. 

By coordinate transformations, other types of wave-fronts (for instance of 
spherical shape) can be studied in a very similar fashion. Outside the support 
of a wave, there is no signal. There the metric space is flat, compatibly with 
the metric tensor when one sets f x = and f y = 0. Basically, for all the 
family of free- waves, we can explicitly provide a full solution to (fT]), including 
situations (not investigated here, but easy to handle) where the polarization 
varies with £. Such a general result extends the analysis developed in [11]. For 
f x and f y constant we get (unbounded) plane waves. Gravitational plane waves 
in vacuum were firstly found in [3] (see also [55], section 35). The results are 
related to a situation different from the one presented here, and the waves look 
actually flat after the introduction of a suitable 'planeness' concept. 

We would like to remark that our analysis has been made possible, in such 
a clean and elegant way, because the electromagnetic tensor U a /3 in @ appears 
with the opposite sign. Our choice can be justified by several reasons. First, 
as we said, with the sign adopted here we can come out with a multitude of 
interesting and significant exact solutions. Actually, by changing the sign of 
the right-hand side in ([42]) . there are a few chances to get reasonable bounded 
solutions, without introducing some passages not well justified from the mathe- 
matical point of view. Secondly, when studying more complex situations where 
the mass tensor in ([2]) plays an effective role, we can interpret Einstein's equation 
as a balance law, between electromagnetic and mechanical energies operating 
in opposition (there might be some analogies with the results in [IB])- This is 
actually what happens for the electromagnetic vortices we are going to study in 
the next sections. 

Old well-known solutions can be also adapted to the new context. This is 
the case for example of the Reissner-Nordstrom metric (see [5_2], section 33.2). 
Relatively to the spherical reference framework (r, 9, </)), we start from the po- 
tential A a — (q/r, 0,0,0). Successively, the nonzero entries of the metric are 
defined as: 

M Q 2 

9oo = 1 - — - — All = -1/900 522 = -r 2 533 = -(rsin6>) 2 (44) 

for some constants M > and Q related to the mass and the charge of the 
black-hole. The difference with the standard case is a switch of a sign in the 
expression of goo (usually this is equal to 1 — M/r + Q 2 /r 2 ). One checks that 
p = and R = (with the exception of the singular point r = 0). Moreover, by 
taking Q = q^/xi^ 2 /2c 4 , the new metric in (j44|) turns out to be solution of the 
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Einstein's equation: 

2 

R af 3 = ^U a0 (45) 

that corresponds to ([TJ) with the right-hand side tensor |(3J), in contrast to the 
classical approach where R a p — — {xp 2 V 'c^Wap (see for instance [3], p. 117). 

It is important to remark that, with such a new setting, one can remove the 
constraint M > 2Q, granting in this way the existence of the horizons at any 
regime, even with small masses. A similar modification of the coefficients (that 
is: q 2 — > — q 2 ) also works in the case of the Kerr-Newman metric (see [H]). 
This also shows that switching the sign of U a p is not a traumatic choice, but 
can give rise to meaningful aspects never investigated before. The possibility of 
dealing with electric objects with small masses reopens the path to the study of 
gravitational electron models, without the use of involved geometrical extensions 
(see [5]). We examine this possibility in section 6. 

Wc finally provide a definition of mass density p m , which is derived from the 
coefficient Moo of the mass tensor in ([6]). In the case when g a p = {1,-1,— 1,-1} 
and | Vb| = c, we set: 

e M o e / p \ , , 

Pm = — = — IP" T J ( 46 ) 

/J,C p V c / 

If E is the support of the wave, the corresponding mass m is obtained by in- 
tegrating p m over S (we have J„ p m s/—g for a generic metric space). In the 
case of a free- wave we have p = 0. If, in addition the wave has compact support 
E, we have: m = J-^Pm = ( e o/p) J^P = ( £ o/a*) J as E • n = 0. In agreement 
with the Gauss's theorem, we are not in presence of classical charges. Thus, a 
pure photon has no charge and mass, even if there are points inside where p and 
p m can be different from zero. In addition, a photon proceeds straightly at the 
speed of light without dissipation and does not affect the environment (except 
for the region directly touched). The curvature of the space-time is altered at 
its passage, but not in a significant way to influence its motion (according to 
(j3"4")) the structure of the geodesies is only modified in the direction orthogonal 
to that of propagation). There is however a latent possibility to generating 
massive charges, and this happens when the wave is trapped in isolated regions 
of space. We analyze this aspect in the coming sections. 

We finally note that the pressure is a potential. Therefore, our scalar p, 
after dimensional adjustment, could be put in relation with a gravitational type 
potential (see also ©), which is zero for a free single photon and definitely 
different from zero when more waves interact. Thanks to ((8|), p is invariant with 
respect to coordinate changes. In [TT] it is observed that, in wave interactions, 
there is a global electromagnetic energy reduction (with respect to the sum of 
the free constituents), which is compensated by the creation of the potential 
energy associated to p. For rotating photons in a bounded region of space, 
generating a gravitational environment, one may impose Vp = at the border 
(no forces outside the object). In the coming sections, with the help of explicit 
computations, we check exactly what happens in this situation. In the end, our 
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approach will turn out to be a reasonable proposal for unifying electromagnetic, 
mechanical and gravitational phenomena (we partly agree with the viewpoints 
expressed in [7J, although our theory is only based on symmetric tensors). 



5 Rotating waves 

We now examine solutions to our model that do not belong to the subspace of 
free-waves. The computations are in general very complicated. Anyway, there 
are some relatively simple configurations where explicit solutions are available 
(see also [TT], chapter 5). 

In this section, the variables are expressed in cylindrical coordinates (r, z, <fi) 
and the solutions will not depend on z. Note that this reference frame is left- 
handed. We first define the potentials. For an integer k > 2 and arbitrary 
constants u > 0, 70 and 71, we set: 

(71 7o' 7i \ 
— Jk+i(ur) sm(cojt — km), , — ; — I Jk+iiuir) cos(cuit — k<f>) 
ui 4fc uj ) 

2 

$ = -^—-^-J k {ujr)cos(cujt-k(j)) (47) 
Zui uj 

for < (j) < 27r, < r < b^jui and any z. In (|4"T)l . by J k , we denote the 
fc-th Bessel function, that is defined to be a solution to the following differential 
equation: 

4{x) + m. _ + Mx) = (48) 

The quantity 5 k turns out to be the first zero of J k , that, for k — 2, takes 
approximately the value (see also table 1): 

8 2 w 5.135622 (49) 

The potentials satisfy the Lorenz gauge condition. From these we get the elec- 
tromagnetic fields: 

E = f-^I. -f kJk{ur) CQS ^ cu) £ _ 0, 7! J' k {ujr) sm(cu>t — k<p) 

\ uj uir 

B = -fo, ^— +liJk{ur)cos(cwt- k<f>), 0) (50) 
c V k I 

In order to check the above expressions it is worthwhile to recall the relations: 
J'k + i(x) + (k + l)Jk+i(x)/x = Jk(x) and J' k (x) - kJ k (x)/x = - J k+ i{x). 

For r = S k /oj, the time-dependent components of E\ and B2 are zero. Note 
also that E and B vanish for r = 0, since J k {x) decays as x k for x — * 0. The 
choice k = 1 is not permitted because it does not allow to prolong with continuity 
the fields up to r = 0. The idea is to simulate a fc-body rotating system in 
equilibrium. This somehow explains why the case k = 1 is not going to produce 
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meaningful solutions. The wave is constrained in a cylinder E, with the magnetic 
field lined up with the axis and the electric field laying on the orthogonal planes. 
One has that divE is constantly equal to 270/w, since the divergence of the time- 
dependent part is zero. Moreover, by defining V = (0, 0, cior/k), one finds out 
that: DV/Dt = (-{cuj/k) 2 r, 0, 0). 

At this point, one can check that all the model equations (f2"2")l . (12"3")) . (12"4")) are 

satisfied, provided we define the pressure as follows: 

P = ^° + — p ~ ~ -2p" J ^ — J fc( wr ) cos («^ - ( 51 ) 

where po is an arbitrary constant. This may be checked by substitution in (1134(1 
after recalling that, using (|48[) with a; = wr, the gradient of p is: 

dp 2c 2 ^ 7o r 2/ry|r/ w 2 r 2 \ 2^7o7i/fc 2 2 \ . . 

g = _2^ Ji(wr)(dn(ajj( _ Jb0) (52) 

Since we have p 7^ 0, we are not in presence of a free-wave. Due to (JSJ) , we 
expect that in the new metric space defined through (p}, the scalar curvature 
R is going to be different from zero, testifying that such a space is effectively 
curved. Indeed, such rotating fronts and bent light rays must be consequences 
of a space-time deformation, where the 4- velocity V a should be compatible with 
the generalized eikonal equation ([7]). Unfortunately, we do not have the metric 
tensor in this complicated situation. In order to get it, one should first create 
from (— $,A) the tensor F a p. Successively, one has to build T a p as in (|2|), 
([3]), ©, where g a p is a generic metric tensor, and finally solve §Q to find the 
unknown g a @. Note that the pressure in (|5ip is equal to po for r = 0, and that 
the stationary components of (E, B, V) form a left-handed triplet. It has to be 
noticed that the momentum of our rotating wave is preserved independently of 
w. high frequencies (w large) are associated with small diameters 5 = Sk/w, and 
vice versa. 

Although we do not have a rigorous proof, the main idea is summarized as 
follows. The k photons rotate around an axis, so that producing, via Einstein's 
equation, a gravitational setting having potential proportional to p (or the scalar 
curvature i?, according to {8])), On the other hand, the metric space where such 
a wave is naturally embedded is such that the light rays follow circular orbits 
(geodesies) and the wave- fronts satisfy an eikonal type equation (see (J7|)). The 
setting may help to better understand the physics of relativistic rotating bodies 
(see for instance [24 ). By breaking the equilibrium, we should get k independent 
free-photons travelling along straight-lines in an almost flat space with p = 
(see section 4). In this passage, Einstein's equation guarantees global energy 
and momentum preservation (see also section 6). 
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Figure 1: Displacement of the left-handed triplet (E,B,V) in the case of a vortex 
ring. Stationary electromagnetic fields are overlapped to oscillating ones, resulting 
in a rotation around the major circumference. The first ones provide for charge and 
magnetic momentum, the second ones are responsible for the curvature of the space 
and the stability of the configuration. 

The cylinder £ is unbounded, so that the energy of our wave turns out 
to be infinite. We can however build more physical solutions, by imposing 
that E is a toroid with major diameter rj and minor diameter S (see figure 
1). The stability of vortex rings, from the fluid dynamical point of view, is 
impressive. We expect the same to happen to our solutions, although, due to 
the difficulty of the problem, such an analysis looks very hard. The numerical 
results in [5] , show that it is possible to build electromagnetic waves constrained 
in toroid structures, ranging from the case of a classical ring to a spherical Hill's 
type vortex, where the inner hole of the ring is reduced to a segment. The 
investigation could help us to understand some electric phenomena, such as ball 
lightning or other self-organized plasmas (see for instance [5]). For the usual 
ring, the study in 8 indicates that, for 77 slightly bigger than 6, the section 
fi of the toroid is practically a circle. Thus, one may consider solutions where 
the electric field belongs to the vertical section of the ring, and the magnetic 
field circulates around the axis (see figure 1). Hence, in good approximation, 
the expression of the fields can be taken equal to that just obtained for the 
case of the cylinder. With this setting, we carry out some computations in the 
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coming section, with the aim to illustrate quantitatively that electrons could be 
modelled by electromagnetic vortex rings of the right size, charge and mass. 



6 A model for the electron 

The idea that matter at atomic level has a structure recalling that of vortex rings 
dates back to [35] • At that time, however, the knowledge of atomic physics was 
very poor. Based on the Maxwell's model, a first effort to describe a spherical 
electron with the help of electromagnetic fields was developed in [23] (see also 
[17]). The project could not be realized due to the scarce information available, 
both theoretical and experimental. For instance, the concept of spin was not 
yet introduced. Moreover, general relativity was developed only ten years later, 
and this, to our opinion, must be a fundamental ingredient for the construction 
and the stability of the model. Successively, a purely geometrical approach was 
considered in [26], with the introduction of the so called geons, that is toroid 
structures in a suitable space-time, roughly approximated by black holes through 
the Kerr-Newman metric (for more recent results, see [5]). The literature has 
many other scattered contributions on the subject. In the end, all the proposed 
models are a bit rough and do not take into account all the features of a particle. 
Nevertheless, since now we exactly know what a photon is, additional elements 
are available for a deeper investigation, and the solution we are going to show 
will be realistic in all the aspects. 

Let us consider figure 1. The volume of E is 2ir 2 rj5 2 , therefore, the total 
charge is given by the integral of p — divE over E: 

f 2 70 47T 2 6 ^ 2 

q = eo / = 7o (53) 

is w w 

Thus, we deduce: 

= quj3 CU\ 
70 47r 2 e ^ 2 47r 2 e ^ 1 ' 

where 5 = 5k /u>, 5k being the first zero of Jk- Successively, by integrating the 
pressure on E, we obtain: 



p = 2irrj 



f c 2 ujj r 2 M7o r f-, uj z r 2 \ 



r drdcf> (55) 



where fi is the circle of radius 5 = 5k/ u>, and we used the stationary part of 
the expression in (|51|) (note that the dynamical part has zero average in fl). 
Proceeding with the computation we get: 



p = 4n 2 ri 



Por 2 c 2 ^7 r 4 
2 4fc 2 



2n 2 np 5l tt 2 ?^^ / c 2 



1 4 

M7o r 


^7 2 r 6 " 


4w 2 


12fc 2 _ 


M7o 




u 3 


3w 3 
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= ~^ + — + le^^T" 1 ) (56) 

where we used (|54|) and we set Ak — S^/k 2 . For convenience, we report in table 
1 the values of 8k and Ak, for various k. 



k 


4 


4/fc 


A k = 8l/k* 


2 


5.1356 


2.567 


6.593 


3 


6.3802 


2.126 


4.522 


4 


7.5883 


1.897 


3.598 


5 


8.7715 


1.754 


3.077 


6 


9.9361 


1.656 


2.742 


7 


11.086 


1.583 


2.508 


8 


12.225 


1.528 


2.335 


9 


13.354 


1.483 


2.201 


10 


14.476 


1.447 


2.095 


11 


15.590 


1.417 


2.008 


12 


16.698 


1.391 


1.936 


13 


17.801 


1.369 


1.875 


14 


18.900 


1.350 


1.822 



Table 1: Values of 5k and associated quantities for various k. 



After recalling the definition (|46p. we now require that p m = for r = 
(zero mass density at the center of the body, in part justified by the fact that p 
is quadratical with respect to r near the origin) . This implies that the constant 
Pq has to be equal to: c 2 p — 2^/qc 2 /uj. We also remark that a pure rotating wave 
with 70 = and 71 7^ (see (|50]>) turns out to have m = / s p m = (since p = 
and p = 0), according to the fact that a photon is uncharged and massless, even 
if it is constrained in a bounded region of space. 

Finally, concerning the mass of our toroid, we can write: 

eo f ( P\ qA k g 2 ( Ak \ , 

m ^ thy?-*) = ~ -ipT ~ t^w It~ V (57) 

Let us assume that 2 < k < 5, so that Ak > 3 (see table 1). If q is positive, the 
above expression is certainly negative and cannot represent a mass. If instead 
q = —e is the electron charge, we have chances to get something interesting. In 
such a situation, the last expression takes the form: 

in = ~7~ ~~ — (^-l)^- (58) 

where a = e 2 /2hceo is the fine structure constant and h is the Planck constant. 
It is important to observe that m does not depend on w and 5 (we discuss this 
aspect later). 
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Successively, we may try to find the major diameter 77. As an extra condition, 
we impose that the radial component of Vp is zero at the boundary of f2. In fact, 
assuming that p is a kind of gravitational potential (see the end of section 4), 
this means that no other external forces are pushing on the boundary (except 
for those related to the centripetal acceleration and the unbalanced Lorentz 
condition E + VxB^O), producing a sort of equilibrium at the surface of E. 
Thus, in order to determine 77, we set to zero the first relation in (|52|) . at the 
point Sk/co. Considering that Jk{$k) = 0, we get the relation: 



Substituting in (|59[) the expression of 70 given in l]54p and introducing the 
constant a, one obtains: 



fj, ah 



e 2ir 2 c \ J±k ' 

We can go back to (f58|) . With the new value of 77, we can finally write: 

' A * (61) 



/J 6(i4 fc - 1) 

Let us discuss the case k = 2. Assuming to know the values of e and to, from 
(|6"l"T) . we deduce the estimate: 

u = 1.295 — w 2.27 x 10 11 Coulomb/Kg (62) 

TO 

A similar estimate of the constant fi was recovered in [TT] by different arguments. 
Successively, from (l60l) we can compute 77: 

127T^TOC 

that agrees quite well with the Compton radius of an electron. 

Let us see what happens to the electromagnetic energy. We take the sta- 
tionary parts of the fields in (f50|) . With computations similar to those carried 
out above, we can evaluate the energy: 



/ (|E| 2 + c 2 |B| 2 ) 



7T£o?7 



r 2 \2 



For k = 2, using the values already computed, we arrive at 

1 / 3 
2-( 1+ 2l 



£ = Ul + 7^7-)mc 2 w .614 toc 2 (65) 
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The remaining part of the energy may be provided by the dynamical part. At 
this point, one computes 71 in (|50p . in order to get a total energy corresponding 
to mc 2 = (eo/fJ,) / s (pc 2 —p), that, written in this form, takes the meaning of a 
potential energy (recall that p is constant in E). 

By analyzing 500 m (|44[) . we can get further information. For example, if 
the radius r] of the particle coincides with the horizon, so that (Q/i]) 2 « 1, one 
should have: x ~ 32(7rc 2 7/e /pe) 2 w 0.145. 

In the end, we were able to prove the following facts. It is possible to build 
electromagnetic vortex rings having the size (2rj), the charge (— e) and the mass 
(m) of an electron (see also [5]). Note that some computations could have been 
made in a different way. For instance, in the definition of p m (see (|46[1 ), we can 
include a further multiplicative adimensional constant. The aim was however 
to show that there exists at least an admissible setting. 

Since our particles are made of rotating photons, it is natural to associate 
an intrinsic spin to them, even if they are not mechanically rotating around 
an axis. There is also a magnetic field inside the toroid, having zero dipole 
moment. This implies that the particle is 'neutral' from the magnetic viewpoint, 
but, under external solicitations, it may exhibit a magnetic momentum. The 
triplet (E, B,V) is left-handed and this has something to do with the concept 
of chirality. In fact, the specular image of our particle, obtained by switching 
the sign of the charge, produces an anti-particle (the positron), sharing the 
same properties of the electron, as it can be proven with similar computations 
(compare (l20l) and (j3Tj) b When an electron and a positron collide, they destroy 
their stable geometrical environment, producing pure photons departing from 
the impact site according to standard mechanical rules. 

There is no a single solution, but a multitude of them, depending on the 
parameter to. As a matter of fact, maintaining the amplitude of 77, we can con- 
struct thin rings £ having the sections fl rotating at very fast speed. Their 
charge and mass densities are high and the volume small. Conversely, one can 
build fat rings with lower angular velocity. Probably, there exists an intermedi- 
ate solution optimizing the energy, but, with an analysis at the first order, we 
were not able to detect it. Thus, we can associate to our particle infinite fre- 
quencies. However, the whole object has a fixed global diameter 2rj, that can be 
formally taken as the wave-length of a suitable fictitious internal clock (see for 
instance [15] ). These facts are not disturbing however. They make the electron 
a very adaptable particle when interacting with the outside world, explaining 
why it can emit photons with a continuous spectrum and escape to finer mea- 
surements. In order to justify these assertions, one has to assume that the bare 
particle is 'dressed up' with a cloud of photons, bringing different frequencies 
and responsible for the transmission of the electric charge to the environment. 
At this point, the discussion becomes more delicate and the interested reader 
is addressed to [TT], chapter 6, and [15], where the author gives his personal 
interpretation on issues regarding the quantization of matter. 

Switching the sign of both the electric and magnetic fields, should bring 
to a positive left-handed particle. Its structure is going to be quantitatively 
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different from that of the electron, since the forces in (j2~4"|) find their balance 
in an alternative way. We do not think however that the standard ring will 
be stable in this case, while we guess that a Hill's type vortex would be more 
appropriate (see [Hj)- On the other hand, if also protons can be described by 
our model, we do expect they assume complicated shapes (see and display 
peculiar substructures, as it is revealed by experiments. 
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